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Abstract 

We calculate the cross section for e + e~ annihilation into three jets for massive 
quarks at next-to-leading order in perturbative QCD, both on and off the Z 
resonance. Our computation allows the implementation of any jet clustering 
algorithm. We give results for the three-jet cross section involving b quarks 
for the JADE and Durham algorithm at cm. energies y/s = mz- We also 
discuss a three-jet observable that is sensitive to the mass of the b quark. 
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Jets of hadrons, which originate from the production and subsequent fragmentation of 
quarks and gluons in high energy electron positron annihilation have been among the key 
predictions [IJQ of quantum chromodynamics (QCD). For precision tests of QCD the e + e~ 
experiments provide a particularity clean set-up. There exist a number of jet observables 
that are well-defined (i.e., infrared-finite) in QCD, and which can be calculated perturba- 
tively as an expansion in the strong coupling a s . The next-to-leading order (NLO) QCD 
corrections to three-jet production were computed more than a decade ago for massless 
quarks, and subsequent implementations |5|-|§ of these results have been widely used for 
tests of QCD with jet physics. 

To date huge samples of jet events produced at the Z resonance have been collected both 
at LEP and SLC. From these data large numbers of jet events involving b quarks can be 
isolated with high purity using vertex detectors. For detailed investigations of b jets quark 
mass effects must be taken into account in the theoretical predictions |TD[. Specifically, 
knowledge of the NLO three-jet fraction for non-zero quark mass opens the possibility to 
measure the mass of the b quark from b jet data at the Z peak JTT|. Further applications 
include precision tests of the asymptotic freedom property of QCD by means of three-jet 



rates measured at various center-of-mass energies, also far below the Z resonance JT2 . 

As far as massive quarks are concerned the three, four, and five jet rates are known 
to leading order (LO) in a s only [pL3|— pT5 1 . In this Letter we report the calculation of the 



e + e annihilation cross section into three jets involving a massive quark antiquark pair at 
next-to- leading order QCD ||IE|,|r?fl. The determination of this cross section cr% LO to order 



'NLO 

a 2 s consists of two parts: First, the computation of the amplitude of the partonic reaction 
e + e~ — > 7*, Z* — > QQg at leading and next-to- leading order in the QCD coupling. Here Q 
denotes a massive quark and g a gluon. We have calculated the complete decay amplitude 
and decay distribution structure for this reaction. This allows for predictions including 
oriented three-jet events. The differential cross section involves the so-called hadronic tensor 



which contains five parity-conserving and four parity-violating Lorentz structures. Second, 
the leading order matrix elements of the four-parton production processes e + e~ — > Z*, 7* — > 
ggQQ, QQqq, QQQQ are needed. Here q denote light quarks which are taken to be massless. 

The calculation of the 0(a 2 ) virtual corrections to the process e + e~ — > QQg with 
massive quarks is straightforward albeit tedious. Non-neglection of the quark mass leads 
to a considerable complication of the algebra. The infrared (IR) and ultraviolet (UV) 
singularities, which are encountered in the computation of the one-loop integrals, are treated 
within the framework of dimensional regularization in D = 4 — 2e space-time dimensions. 
We remove the UV singularities by the standard MS renormalization. An essential aspect of 
our computation is to show that the IR singularities of the virtual corrections are cancelled 
by the singularities resulting from phase space integration of the squared tree amplitudes 
for the production of four partons. Different methods to perform this cancellation have 
been developed (see []7|-|j| and references therein). We use the so-called phase space slicing 
method elaborated in @ . The basic idea is to "slice" the phase space of the four parton final 
state by introducing an unphysical parton resolution parameter s min <C sy cut , where y cut is 
the jet resolution parameter. The parameter s min splits the phase space into a region where 
all four partons are "resolved" and a region where at least one parton remains unresolved. 
For massless partons, the resolved region may be conveniently defined by the requirement 
that all invariants Sy = (fcj + kj) 2 built from the parton momenta fc, are larger than the 



parameter s m i n . We have modified this definition slightly to account for masses [ |I~6|j . 

In the unresolved region soft and collinear divergences reside, which have to be isolated 
explicitly to cancel the singularities of the virtual corrections. This is considerably simplified 
due to collinear and soft factorizations of the matrix elements which hold in the limit 
Smin — > 0. In the presence of massive quarks, the structure of collinear and soft poles is 
completely different as compared to the massless case. As an example, we would like to 
discuss for e + e~ — > Q(k 1 )Q(k 2 )g(k 3 )g(k4) the limit where one gluon, say g(k^), becomes soft. 



In this limit, the squared matrix element can be written as a universal factor multiplying 
the squared Born matrix element for e + e~ — > QQg. The integration over the soft gluon 
momentum fc 4 can then be carried out analytically in D dimensions. In the soft limit k± — > 
the squared matrix element reads 
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where g s is the strong coupling constant, Nq 
m denotes the quark mass. Each of the three terms in (p]) can now be integrated over the 
appropriate soft phase space volume which we define by the conditions t a ^ + t u < 2 s m i n 
(leading color terms), and £14 +£24 < 2 s m i n (subleading color term). The complete soft fac- 
tor S(ki, k 2 , k 3 ) multiplying the squared Born matrix element for e + e~ — > Q(k 1 )Q(k 2 )g(k 3 ) 
which is obtained by this integration reads: 
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where (3 — Jl — Am 2 / (£ 12 + 2m 2 ), lu = (1— /3)/(l+/3), and /1 is an arbitrary scale introduced 
to keep a s dimensionless in D dimensions. The poles in e exhibited in (0) (and additional 
poles from the collinear region of phase space which we do not show explicitly) can now be 
cancelled against the IR poles of the one-loop integrals entering the virtual corrections. One 
is then left with a completely regular differential three-parton cross section which depends 



on s r 



The contribution to cr% LO of the "resolved" part of the four-parton cross section is finite 
and may be evaluated in D = 4 dimensions, which is of great practical importance. It also 
depends on s min and is most conveniently obtained by a numerical integration. Since the 
parameter s mi - n is completely arbitrary, the sum of all contributions to cr% LO (y cut ) must not 
depend on s m i n . In the soft and collinear approximations one neglects terms which vanish 
as s m i n 0. This limit can be carried out numerically. Since the individual contributions 
depend logarithmically on s m j n , it is a nontrivial test of the calculation to demonstrate that 
a NLO becomes independent of s min for small values of this parameter. Moreover, in order 
to avoid large numerical cancellations, one should determine the largest value of s m i n which 
has this property. 

The three jet cross section depends on the experimental jet definition. We consider here 



the JADE [0 and Durham |19| clustering algorithms, although other schemes [|20|] can also 
be easily implemented. We have checked that we recover the result of || in the massless 
limit. 

We now present our results for the cross section cr^ b LO for b quarks in the MS scheme. 
We require that at least two of the jets that remain after the clustering procedure contain 
a b or b quark f21 ]. We use the b quark mass parameter rrib(fi) defined in the MS scheme 



at a scale /i. The asymptotic freedom property of QCD predicts that this mass parameter 



decreases when being evaluated at a higher scale. With mb(rrib) = 4.36 GeV |22j and a s {mz) 
= 0.118 |23[] as an input and employing the standard renormalization group evolution of the 
coupling and the quark masses, we use the value mf,(/i = mz) = 3 GeV. 

Figs. la,b show the three jet cross section a^ b LO at s/s = \i = mz with b quarks of mass 
ruf, = 3 GeV as a function of y m in — s m i n / {sy C ut) for the JADE and Durham algorithms at 
a value of the jet resolution parameter y cut = 0.03. It can be clearly seen that the cross 
section reaches a plateau for small values of the parameter y m in- The error of the numerical 
integration becomes bigger as y m in — » 0. In order to keep this error as small as possible 



without introducing a systematic error from using the soft and collinear approximations, we 
take y m i n = 10~ 2 for the JADE algorithm and y m in = 5 x 10~ 3 for the Durham algorithm. 
These values are used in Figs. lc,d, where we plot ctjvlo as a function of y cut together 
with the LO result. The QCD corrections to the LO result are quite sizable, as known also 
in the massless case. The renormalization scale dependence which is also shown in Figs. 
lc,d is modest in the whole y cu t range exhibited for the Durham and above y cu t ~ 0.01 for 
the JADE algorithm. Below this value perturbation theory is not applicable in the JADE 
scheme. 

The effects of the b quark mass at the Z peak may be exhibited with the following ratio 

m 



{Vcut) ~ Rf sc (y C u t y { ) 

Here we define -R3 = cr 3,b /a(e + e~ — > bb), and R^ dsc is the three-jet fraction for the four light 
quarks with no flavor tagging. (Note that B 7^ 1 in the limit — > due to the different 
definitions of and R^ dsc .) The LO and NLO results for the observable B are shown in 
Figs. 2a,b for the JADE and the Durham algorithms. We took the massless 0(a 2 s ) results 
from ||20|. As both the LO and NLO terms in R\ depend on nib it is clear that comparison 



of our result with measured values of B{y cut ) would allow for an unambiguous determination 
of mb within a given renormalization scheme. In view that the b mass effect in B - and 



in other observables, for instance the differential two-jet distribution [16j - is only of the 
order of a few percent this constitutes an experimental challenge. Moreover, further work 
is needed to assess in detail the theoretical uncertainties involved [^J . Yet such an analysis 
would be worth the effort: it would be the first determination of the b quark mass at a high 
scale, and it might also experimentally establish the "running" of a quark mass as predicted 
by QCD. 

Summarizing we have computed the NLO QCD corrections for e + e~ — > 3 jets for mas- 



sive quarks. Our results, which we shall report on in detail in future work, should find 
applications to a number of precision tests of QCD involving b and c quarks at various cm. 
energies, and to theoretical investigations of top quark production at very high-energetic 
electron positron collisions. 
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Figure 1: Figs, la and lb show ct^lo as defined in the text at y/s = /x = mz as a function 
of Vmin = s m i n /{sy cu t) for the JADE and Durham algorithms at a value of the jet resolution 
parameter y cut = 0.03 with m&(/i = mz) = 3 GeV and a s (/x = mz) = 0.118. Figs, lc and 
Id show a 3 ' b as a function of y cu t for the JADE and Durham algorithms, respectively. The 
dashed line is the LO result. The NLO results are for ji = mz (solid line), \x = (dotted 
line), and fi = 2m z (dash-dotted line). Initial-state photon radiation is not included in the 
cross sections. 
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a) JADE algorithm 
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b) DURHAM algorithm 
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Figure 2: The ratio B of eq. (||) as a function of y cu t for (a) the JADE and (b) the Durham 
algorithm, using values for nib and a s as in Fig. 1. The dashed line is the LO result. The 
points and the solid line are the NLO result for ix = mz- For illustrative purposes, the 
dash-dotted line shows the LO result for = 5 GeV. 
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